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Abstract 

Inclusive jet production is studied in neutral current deep-inelastic positron-proton scat- 
tering at large four momentum transfer squared Q 2 > 150 GeV 2 with the HI detector at 
HERA. Single and double differential inclusive jet cross sections are measured as a func- 
tion of Q 2 and of the transverse energy Et of the jets in the Breit frame. The measure- 
ments are found to be well described by calculations at next-to-leading order in perturba- 
tive QCD. The running of the strong coupling is demonstrated and the value of a s {Mz) 
is determined. The ratio of the inclusive jet cross section to the inclusive neutral cur- 
rent cross section is also measured and used to extract a precise value for ot s (Mz) = 
0.1193 ± 0.0014 (exp.) t° Q Q f Q (th.) ± 0.0016 (pdf) . 
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1 Introduction 



Jet production in neutral current (NC) deep-inelastic scattering (DIS) at HERA provides an 
important testing ground for Quantum Chromodynamics (QCD). The Born contribution in DIS 
(figure [T^) gives only indirect information on the strong coupling a s via scaling violations of the 
proton structure functions. At leading order (LO) in a s additional processes contribute: QCD- 
Compton (figure [TJ5) and boson-gluon fusion (figure [Q;). In the Breit frame of reference [1,2], 
where the virtual boson and the proton collide head on, the Born contribution generates no 
transverse momenta. Partons with transverse momenta are produced in lowest order by the 
QCD-Compton and boson-gluon fusion processes. Jet production in the Breit frame therefore 
provides direct sensitivity to a s and allows for a precision test of QCD. 

Analyses of inclusive jet production in DIS at high four momentum transfer squared Q 2 
were previously performed by the HI [2] and ZEUS [3,4] collaborations at HERA. Perturbative 
QCD (pQCD) calculations supplemented with hadronisation corrections were found to describe 
the data. The strong coupling a s and the gluon density in the proton were both extracted. 

In this paper new measurements of the inclusive jet cross section are presented, based on 
data corresponding to twice the integrated luminosity and a higher centre-of-mass energy than 
in the previous HI analysis [2]. The larger data set together with improved understanding 
of the hadronic energy measurement significantly reduces the total uncertainty of the results. 
Differential inclusive jet cross sections are measured as functions of the hard scales Q 2 and 
the transverse jet energy E T in the Breit frame in the ranges 150 < Q 2 < 15000 GeV 2 and 
7 < E T < 50 GeV. In addition, the ratio of the jet cross section to the inclusive NC DIS cross 
section, in the following referred to as the normalised inclusive jet cross section, is determined. 
This observable benefits from a partial cancellation of experimental and theoretical uncertain- 
ties. The measurements are compared with pQCD predictions at next-to-leading order (NLO), 
and the strong coupling a s is determined from a fit of the predictions to the data. 

2 Experimental Method 

The data were collected with the HI detector at HERA in the years 1999 and 2000. Dur- 
ing this period HERA collided positrons of energy E e = 27.5 GeV with protons of energy 



Figure 1 : Deep-inelastic lepton-proton scattering at different orders in a s : (a) Born contribution 
C(l), (b) QCD Compton scattering 0(a s ) and (c) boson-gluon fusion 0{a s ). 
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E p = 920 GeV giving a centre-of-mass energy y/s = 319 GeV. The data sample used in this 
analysis corresponds to an integrated luminosity of 65.4 pb -1 . 

2.1 HI detector 

A detailed description of the HI detector can be found in [5, 6]. HI uses a right-handed co- 
ordinate system with the origin at the nominal interaction point and the z- axis along the beam 
direction, the +z or "forward" direction being that of the outgoing proton beam. Polar angles 
9 and azimuthal angles <fi are defined with respect to this axis. The pseudorapidity is related to 
the polar angle 9 by 77 = — lntan(6 l /2). The detector components important for this analysis 
are described below. 

The electromagnetic and hadronic energies are measured using the Liquid Argon (LAr) ca- 
lorimeter in the polar angular range 4° < 9 < 154° and with full azimuthal coverage. The 
LAr calorimeter consists of an electromagnetic section (20 to 30 radiation lengths) with lead 
absorbers and a hadronic section with steel absorbers. The total depth of both sections varies 
between 4.5 and 8 interaction lengths. The energy resolution is <Je/E = 12%/ \J~E~j GeV ©1% 
for electrons and oeIE = 50%/ \f E / GeV © 2% for hadrons, as obtained from test beam 
measurements [7]. In the backward region (153° < 9 < 177°) energy is measured by a 
lead/scintillating fibre Spaghetti-type Calorimeter (SpaCal) composed of an electromagnetic 
and a hadronic section. The energy resolution of the SpaCal is a E /E « 7% / \J~E / GeV © 1% 
for electrons [8]. The central tracking system (20° < 9 < 160°) is located inside the LAr 
calorimeter and consists of drift and proportional chambers, complemented by a silicon vertex 
detector [9] covering the range 30° < 9 < 150°. The chambers and calorimeters are surrounded 
by a superconducting solenoid providing a uniform field of 1.16 T inside the tracking volume. 

The scattered positron is identified as an electromagnetic cluster in the LAr calorimeter 
with an associated track. The remaining clusters in the calorimeters and charged tracks are 
attributed to the hadronic final state which is reconstructed using an energy flow algorithm 
that avoids double counting of energy. The luminosity is determined by measuring the Bethe- 
Heitler process (ep — > epj), where the photon is detected in a calorimeter close to the beam 
pipe at z = —103 m. 

2.2 Event and jet selection 

NC DIS events are selected by requiring the scattered positron to be detected in the LAr 
calorimeter with a reconstructed energy E' e exceeding 11 GeV and a polar angle 9 C < 153°. 
These requirements ensure a trigger efficiency of greater than 98%. The ^-coordinate of the 
event vertex is required to be within ±35 cm of the average position of the interaction point. 
This condition reduces contributions from beam induced background and cosmic muons. Non- 
ep background is further reduced by requiring an event timing which matches the HERA bunch 
crossing. The total longitudinal energy balance must satisfy 45 < ~ Pz,i) < 65 GeV, 

where the sum runs over all detected particles. This requirement reduces the contributions of 
the photoproduction background and of DIS with initial state photon radiation for which the 
escaped positron or photon in the — ^-direction leads to values of this observable lower than 
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the expectation 2E e = 55 GeV, for events with losses only along the outgoing proton beam. 
Elastic QED Compton and lepton pair production processes are suppressed by rejecting events 
containing a second isolated electromagnetic deposit and no hadronic activity. The remain- 
ing photoproduction background is estimated using Monte Carlo simulations and found to be 
negligible in all Q 2 and jet Et bins. 

The DIS phase space covered by this analysis is defined by 

150 < Q 2 < 15000 GeV 2 , 
0.2 < y < 0.7 , 

where y, ,quantifies the inelasticity of the interaction. These two variables are reconstructed 
from the four momenta of the scattered positron and the hadronic final state particles using the 
electron-sigma method [10]. 

The jet analysis is performed in the Breit frame. The boost from the laboratory system to the 
Breit frame is determined by Q 2 , y and the azimuthal angle of the scattered positron. Particles 
of the hadronic final state are clustered into jets using the inclusive k T algorithm [11] with the 
p T recombination scheme and with distance parameter R = 1 in the r/-0 plane. The cut — 1.0 < 
r] Lab < 2.5 ensures that jets are well contained within the acceptance of the LAr calorimeter, 
where r/ Lab is the jet pseudorapidity in the laboratory frame. Every jet with 7 < Et < 50 GeV 
contributes to the inclusive jet cross section, regardless of the jet multiplicity in the event. In 
total 23714 jets pass the analysis cuts. 

In addition, the normalised inclusive jet cross section is investigated, calculated as the ratio 
of the number of jets to the number of selected NC DIS events in the y range defined above. 
This observable equals the average jet multiplicity of NC DIS events within the given phase 
space. Jet cross sections and normalised jet cross sections are studied as a function of Q 2 and 
Et- 

2.3 Cross section determination 

In order to extract the cross sections at hadron level, the experimental data are corrected for 
limited detector acceptance and resolution. The correction factors are determined using simu- 
lated NC DIS events. The generated events are passed through a detailed simulation of the HI 
detector and subjected to the same reconstruction and analysis chain as the data. The following 
Monte Carlo event generators are used for this purpose: DJANGOH [12] using the Color Dipole 
Model as implemented in ARIADNE [13], and RAPGAP [14] using matrix elements matched 
with parton showers. Both RAPGAP and DJANGOH provide a good description of the data 
in both the inclusive and the jet sample. The purity of the jet sample, defined as the fraction 
of events reconstructed in a bin that originate from that bin on hadron level, is found to be 
larger than 60% in all analysis bins. Correction factors are determined as the ratio of the cross 
section obtained from particles at hadron level to the cross section calculated using particles 
reconstructed in the detector. This correction is applied bin-by-bin in Q 2 and Et- Arithmetic 
means of the correction factors determined by RAPGAP and DJANGOH are used, and half of 
the difference is assigned as model uncertainty. The correction factors deviate typically by less 
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than 20% from unity. The effects of QED radiation are corrected for using the HERACLES [15] 
program. The size of these corrections is typically 10% for the jet cross sections and 5% for 
the normalised jet cross sections. 

For the normalised jet cross sections the ratio of the number of jets to the number of NC 
DIS events is calculated on detector level, and it is this ratio which is corrected for detector and 
QED effects. 

2.4 Systematic errors 

The following sources of systematic uncertainty are considered: 

• The positron energy uncertainty is 0.7% to 3% depending on the ^-impact point of the 
positron in the calorimeter. Uncertainties in the positron reconstruction affect the event 
kinematics and thus the boost to the Breit frame. The resulting uncertainty on the cross 
sections and normalised cross sections is typically 0.5%. 

• The positron polar angle systematic uncertainty is between 1 and 3 mrad. The resulting 
uncertainty on the cross sections and normalised cross sections is typically 0.5%. 

• The energy scale uncertainty of the reconstructed hadronic final state is estimated to be 
2%, dominated by the uncertainty of the LAr hadronic energy scale. The resulting un- 
certainty on the cross sections and normalised cross sections is typically in the range 
1 to 4%. 

• The luminosity measurement uncertainty leads to an overall normalisation error of 1.5% 
for the jet cross sections. 

• The model dependence of the data correction is estimated as described in section [2731 It 
is below 10% in most of the bins and typically 2%. 

• An error of 1% is estimated from the uncertainty of the QED radiative correction [16]. 

The dominant experimental uncertainties on the jet cross section arise from the model depen- 
dence of the data correction and from the LAr hadronic energy scale uncertainty. The individual 
contributions are added in quadrature to obtain the total systematic uncertainty. The correlations 
of the errors among the different bins are treated using the same procedure as described in [2]. 
The uncertainties of the luminosity measurement and of the positron polar angle are each as- 
sumed to be fully correlated between the bins. The error on the positron polar angle and the 
QED radiative corrections is assumed to be uncorrelated. The remaining sources of system- 
atics, namely the positron energy scale, the hadronic final state energy scale and the model 
dependence are equally shared between correlated and uncorrelated parts. For the normalised 
jet cross sections systematic uncertainties are reduced and the luminosity uncertainty cancels. 
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3 NLO QCD Calculation 



Reliable quantitative predictions of jet cross sections in DIS require the perturbative calculations 
to be performed to at least next-to-leading order of the strong coupling. In order to compare with 
data, hadronisation corrections have to be applied to the perturbative calculations. By using the 
inclusive kx jet algorithm [11] the observables in the present analysis are infrared and collinear 
safe and the hadronisation corrections are small. In addition, by applying this algorithm in the 
Breit frame, jet cross sections can be calculated in pQCD, since initial state singularities can be 
absorbed in the definition of the proton parton densities. 

The theoretical prediction for the jet cross section is obtained using the NLOJET++ pro- 
gram [17], which performs the matrix element integration at NLO of the strong coupling, 
0(a 2 ). The strong coupling is taken as a s (Mz) = 0.118 and is evolved as a function of the 
renormalisation scale at two loop precision. The calculations are performed in the MS scheme 
for five massless quark flavours. The parton density functions (PDFs) of the proton are taken 
from the CTEQ6.5M set [18]. The factorisation scale /if is chosen to be Q and the renormalisa- 
tion scale /i r is chosen to be the E T of each jet. Running of the electromagnetic coupling with 
Q 2 is taken into account. No QED radiation is included in the calculation since the data are 
corrected for this effect. Electroweak effects due to Z° boson exchange are determined using 
the LEPTO event generator [19] and are applied as correction factors to the calculation. 

The hadronisation correction factor is calculated for each bin as the ratio of the cross section 
defined at hadron level to the cross section defined at parton level. These correction factors are 
determined with the same Monte Carlo event samples used to correct the data from detector 
to hadron level. The correction factors applied to the perturbative calculations are calculated 
as the average of the values from DJANGOH and RAPGAP, as described in section 12.31 The 
hadronisation correction factors differ typically by less than 10% from unity and agree at the 
level of 2% between the two Monte Carlo simulations. 

The theory uncertainty includes the hadronisation correction error and the uncertainty re- 
lated to the neglected higher orders in the perturbative calculation. The systematic error at- 
tributed to the hadronisation correction is taken to be half of the difference between the correc- 
tion factors obtained using RAPGAP and DJANGOH. The dominant uncertainty is related to 
the NLO accuracy and is estimated by a variation of the chosen scales for ji r and fif by arbitrary 
but conventional factors in the range from 0.5 to 2 applied to the nominal scales. In seven out 
of the 24 bins in Q 2 and E T the dependence of the pQCD calculation on fi r is not monotone, 
i.e. the largest deviation from the central value is found for factors within the range 0.5 to 2. In 
such cases the difference between maximum and minimum cross sections found in the variation 
interval is taken, in order not to underestimate the scale dependence. Over the whole phase 
space, the uncertainty due to the renormalisation scale is found to be at least a factor of three 
larger than that due to the factorisation scale. The contributions from both scale variations are 
added in quadrature. 

In order to calculate the normalised inclusive jet cross sections, the prediction of the in- 
clusive jet cross section is divided by the prediction of the NC DIS cross section. The latter 
is calculated at NLO, 0(a s ), with the DISENT package [20], using the same settings as for 
NLOJET++ and with the renormalisation and factorisation scales set to Q. Again, the scale 



8 



uncertainties are determined by independent variations of /i r and fif in the range from 0.5 to 2 
around the nominal value. The scale uncertainties from the jet and the NC DIS part are assumed 
to be uncorrelated. Consequently, the scale uncertainty for the ratio is estimated by adding both 
contributions in quadrature. If the uncertainties are assumed to be anti correlated, which leads 
to the largest change, the resulting theory error increases only slightly by a factor of 1.15. The 
uncertainty originating from the PDFs is also taken into account. The CTEQ6.5M set of parton 
densities provides variations which quantify the uncertainty of the central set. The PDF uncer- 
tainties are propagated into the pQCD prediction of the inclusive jet cross section and the NC 
DIS cross section. 

The strong coupling is determined by repeating the perturbative calculations for many values 
of a s (M z ) until the best match of data and theory is found. With NLOJET++ and DISENT 
these calculations are time consuming. A considerable gain in computational speed is provided 
by the fastNLO package [21], which uses a two step strategy to reduce the calculation time. In 
the first step, the integration of the matrix elements is performed, which is the time consuming 
part of the calculation. This step relies for the present analysis on NLOJET++ and DISENT 
and is independent of a s (Mz), PDFs and the renormalisation scale. In the second step, the 
cross sections are calculated with these parameters specified. The interpolations involved in 
this procedure yield a precision of better than 0.2% on the cross section. All theory calculations 
shown in the following are obtained using fastNLO. 

4 Results 

In the following, the differential cross sections are presented for inclusive jet production and 
for normalised inclusive jet production. Tables Q] and [2] list the measured cross sections to- 
gether with their experimental uncertainties, separated into bin-to-bin correlated and uncorre- 
lated parts. These measurements are subsequently used to extract the strong coupling a s . 

4.1 Cross section measurements compared to NLO predictions 

The measured cross sections, corrected for detector and radiative QED effects, are presented 
as single and double differential distributions in figures [2l|4] The data points are shown at the 
average value of the Q 2 or Et in each bin. The results are compared to the perturbative QCD 
predictions in NLO with a s {Mz) = 0.118, taking into account hadronisation effects and Z° 
boson exchange as explained in section [3] 

The single differential inclusive jet cross sections, defined for events with inelasticity 0.2 < 
y < 0.7 and jets with pseudorapidity —1.0 < r] Lah < 2.5, are shown in figure[2]as functions of 
Q 2 and Et- A good description of the data by the theory calculation is observed. 

The double differential inclusive jet cross section is shown in figure [3] as a function of E T 
in six Q 2 bins in the range 150 < Q 2 < 15000 GeV 2 . The data are well described by the 
theory over the full E T and Q 2 ranges, with x 2 /ndf = 16.7/24, taking only experimental errors 
into account. The agreement is also good when Q instead of Et is used in the calculation as 
renormalisation scale (x 2 /ndf = 24.0/24). 
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For NC DIS events in the range 0.2 < y < 0.7 and in a given Q 2 bin the normalised inclusive 
jet cross section is defined as the average number of jets within —1.0 < 7] Lsb < 2.5 per event. 
Figure S] shows the normalised inclusive jet cross section as a function of E T in six Q 2 bins. The 
NLO calculation gives a good description of the data in the full E T and Q 2 range. Compared 
with the inclusive jet cross section, the normalised inclusive jet cross section exhibits a smaller 
experimental uncertainty. 

4.2 Extraction of the strong coupling 

The QCD predictions for jet production depend on a s and on the gluon and the quark density 
functions of the proton. Using the present jet cross section measurements and the parton den- 
sity functions from global analyses of inclusive deep-inelastic scattering and other data, a s is 
determined. 

QCD predictions of the jet cross sections are calculated as a function of a s {jj, r = Et) with 
the fastNLO package. The cross sections are determined using the CTEQ6.5M proton PDFs and 
hadronisation correction factors as described in section[3l Measurements and theory predictions 
are used to calculate a x 2 ( a s) with the Hessian method, where parameters representing system- 
atic shifts of detector related observables are left free in the fit. The experimental shifts (model 
dependence of the correction factors, positron energy scale, positron azimuth, hadronic final 
state energy scale and luminosity) found by the fit are consistent with the quoted uncertainties. 
This method fully takes into account correlations of experimental uncertainties [22] and has 
also been used in global data analyses [23,24] and in previous HI publications [2,25], where a 
detailed description can be found. The experimental uncertainty of a s is defined by that change 
in a s which gives an increase in x 2 of one unit with respect to the minimal value. The theory 
error is estimated by adding in quadrature the deviation of a s from the central value when the 
fit is repeated with independent variations of the renormalisation scale, the factorisation scale 
and the hadronisation correction factor. 

First, individual fits of a s to each of the 24 measurements of the double differential inclusive 
jet cross sections (presented in figure [3]) are made. The resulting o: s (Et) are shown in figure [51 
for all bins. These determinations demonstrate the property of asymptotic freedom of QCD 
and are in agreement with the predicted scale dependence of a s . The a s values at the scale 
E T can also be related to the value of the strong coupling at the Z° mass a s (M z ) using the 
renormalisation group equation at two loops. All 24 measurements are then used in a common 
fit of the strong coupling, which yields 

a s (M z ) = 0.1179 ± 0.0024 (exp.) ±8:0033 ( th ± 0.0028 (pdf) , (1) 

with a fit quality: xV n df = 20.2/23. The dominating experimental uncertainty is due to the 
LAr energy scale and the model dependence of the detector corrections. The renormalisation 
scale variation is the main contribution to the theory uncertainty, which dominates the overall 
uncertainty of this a s determination. The fit is repeated with Q instead of E T as an alternative 
choice of renormalisation scale. It yields a larger but compatible value of the strong coupling 
a s (M z ) = 0.1230 ± 0.0028 (exp.) ±°;gg^ (scale) with x 2 /ndf = 25.2/23. The quoted scale 
error corresponds to the variation of the renormalisation scale as described in section |3] 
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The global fit of the CTEQ6.5M PDFs was made assuming a s (M z ) = 0.118. In order 
to test whether this value of a s (M z ) biases the results obtained using the nominal method 
presented above, a method, similar to the one used in [3], is employed using the PDFs from the 
CTEQ6AB series, which were obtained from global fits assuming different values for a s (Mz). 
The cross section as a function of the strong coupling is interpolated with a polynomial and 
this interpolation is used to determine the best fit of the strong coupling to the data. The result 
obtained with this alternative fit method is found to be compatible within 0.3 standard deviations 
of the experimental error with the value from the nominal method. Hence there is no indication 
for a bias due to the value of the strong coupling assumed for the CTEQ6.5M PDFs. 

The measurements of the normalised inclusive jet cross section are also used to extract the 
strong coupling using the nominal fit method. The resulting ol s {Et) are shown in figure [6l for 
all bins. As the results are consistent over the whole range of Q 2 and E T , combined fits are 
made to groups of data points. To study the scale dependence of a s , the six data points with 
different Q 2 at a given Et are used together, and four values of ol s {Et) are extracted. The 
results are shown in figure [7^, where the running of the strong coupling is also clearly observed. 
Finally, all 24 measurements are used in a common fit of the strong coupling, which yields 

a s (M z ) = 0.1193 ± 0.0014 (exp.) (th.) ± 0.0016 (pdf) , (2) 

with a fit quality of % 2 /ndf = 28.7/23. This result is compatible within errors with the value 
from the inclusive jet cross sections quoted in [Q The normalisation gives rise to cancellations 
of systematic effects, which lead to improved experimental and PDF uncertainties. This de- 
termination of a s (Mz) is consistent with the world average a s [Mz) = 0.1176 ± 0.0020 [26] 
and with the previous HI determination from inclusive jet production measurements [2]. In 
figure the running of the strong coupling is studied using the alternative scale Q instead of 
Et'- the four data points at a given Q 2 are used together, and six values of a s (Q) are extracted. 

The dominating theory error can be reduced at the expense of a larger experimental un- 
certainty by restricting the data points included in the fit to those at higher values of Q 2 . The 
smallest total uncertainty is obtained by a combined fit of the normalised inclusive jet cross 
section for 700 < Q 2 < 5000 GeV 2 , 

a s (M z ) = 0.1171 ± 0.0023 (exp.) t° S (th-) ± 0.0010 (pdf) , (3) 

with a fit quality of x 2 /ndf = !- 2 /3. 



5 Conclusion 

Measurements of inclusive jet cross sections in the Breit frame in deep-inelastic positron-proton 
scattering in the range 150 < Q 2 < 15000 GeV 2 are presented, together with the normalised 
inclusive jet cross sections, defined as the ratio of the inclusive jet cross section to the NC DIS 
cross section within the given phase space. Calculations at NLO QCD, corrected for hadroni- 
sation effects, provide a good description of the single and double differential cross sections as 
functions of the jet transverse energy E T and Q 2 . The strong coupling a s is determined from a 
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fit of the NLO prediction to the measurements. The experimentally most precise determination 
of a s (M z ) is derived from the measurement of the normalised inclusive jet cross section: 

a s (M z ) = 0.1193 ± 0.0014 (exp.) (th.) ± 0.0016 (pdf) . 

Additionally, the PDF uncertainty is significantly reduced compared to the determination from 
the inclusive jet cross section. The dominating source of error is the renormalisation scale 
dependence which is used to estimate the effect of missing higher orders beyond NLO in the 
pQCD prediction. This result shows a level of experimental precision competitive with a s 
determinations from other recent jet production measurements at HERA [27] and those from 
e + e~ data [28] and is in good agreement with the world average. 
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bin number corresponding Q 2 range 



1 


150 < Q 2 


< 200 GeV 2 


bin letter corresponding E T range 


2 


200 < Q 2 


< 270 GeV 2 


a 


7 < E T < 11 GeV 


3 


270 < Q 2 


< 400 GeV 2 


b 


11 < E T < 18 GeV 


4 


400 < Q 2 


< 700 GeV 2 


c 


18 < E T < 30 GeV 


5 


700 < Q 2 


< 5000 GeV 2 


d 


30 < E T < 50 GeV 


6 


5000 < Q 2 


< 15000 GeV 2 







Inclusive jet cross section in bins of Q 2 and E T 



total total single contributions to correlated uncertainty 



bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model dep. 


positron 


positron 


HFS hadr. 


hadronis. 




section 


uncert. 


uncertainty 


uncertainty 


uncertainty 


detector corr. 


energy scale 


polar angle 


energy scale 


correct. 




(in pb) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


factor 


1 a 


73.81 


2.1 


6.8 


4.9 


4.7 


3.5 


0.7 


0.4 


2.6 


1.076 


1 b 


32.44 


3.1 


7.7 


5.8 


5.1 


2.7 


0.8 


0.3 


3.9 


1.035 


1 c 


6.40 


7.0 


10.4 


8.8 


5.6 


1.7 


0.3 


0.6 


5.1 


1.032 


1 d 


0.94 


18.9 


21.0 


19.9 


6.5 


0.8 


1.1 


0.9 


6.1 


1.065 


2 a 


58.06 


2.2 


6.3 


4.6 


4.3 


3.1 


0.1 


0.4 


2.6 


1.075 


2 b 


28.85 


3.1 


9.5 


7.0 


6.5 


4.9 


0.3 


0.1 


3.9 


1.034 


2c 


6.16 


6.8 


10.5 


8.8 


5.7 


1.0 


0.7 


0.0 


5.4 


1.040 


2 d 


0.85 


18.9 


21.7 


20.3 


7.5 


2.2 


1.2 


0.8 


6.9 


1.044 


3 a 


55.16 


2.2 


5.6 


4.1 


3.8 


2.6 


0.2 


0.1 


2.4 


1.085 


3 b 


30.45 


2.9 


8.7 


6.4 


5.9 


4.2 


0.1 


0.1 


3.8 


1.032 


3c 


7.87 


6.0 


10.7 


8.6 


6.3 


3.1 


0.1 


0.1 


5.3 


1.029 


3d 


0.69 


18.1 


21.5 


19.9 


8.2 


2.0 


0.7 


1.1 


7.7 


1.039 


4a 


48.50 


2.3 


5.0 


3.8 


3.3 


1.8 


0.2 


0.5 


2.4 


1.093 


4 b 


26.81 


3.0 


8.0 


5.9 


5.3 


3.8 


0.2 


0.4 


3.5 


1.035 


4c 


8.46 


5.4 


10.7 


8.4 


6.6 


3.1 


0.4 


0.1 


5.6 


1.025 


4d 


1.69 


13.3 


16.4 


14.9 


6.9 


3.2 


0.5 


0.5 


5.9 


1.035 


5 a 


43.02 


2.4 


5.2 


3.9 


3.5 


1.7 


0.2 


0.8 


2.5 


1.103 


5 b 


30.23 


2.9 


6.2 


4.7 


4.0 


2.5 


0.4 


0.7 


2.7 


1.040 


5c 


11.88 


4.5 


13.6 


10.1 


9.1 


7.9 


0.2 


0.6 


4.3 


1.038 


5 d 


2.63 


10.3 


16.5 


13.7 


9.2 


6.3 


0.3 


0.9 


6.5 


1.046 


6 a 


1.79 


10.8 


12.8 


11.8 


5.0 


0.8 


4.0 


0.3 


1.8 


1.083 


6 b 


1.23 


13.4 


22.7 


18.5 


13.1 


10.8 


6.4 


1.9 


2.2 


1.050 


6c 


0.76 


17.8 


27.2 


22.9 


14.7 


10.8 


9.2 


2.0 


2.7 


1.029 


6d 


0.44 


26.8 


34.3 


30.8 


15.2 


14.5 


3.0 


1.0 


3.3 


1.029 



Inclusive jet cross section in bins of E T 



a 


281.43 


0.9 


5.7 


4.0 


4.1 


2.8 


0.2 


0.3 


2.5 


1.084 


b 


150.22 


1.1 


7.7 


5.5 


5.4 


3.8 


0.2 


0.2 


3.6 


1.035 


c 


41.70 


2.1 


9.6 


7.0 


6.6 


4.1 


0.1 


0.2 


5.1 


1.033 


d 


7.29 


5.1 


12.1 


9.4 


7.6 


3.8 


0.5 


0.3 


6.5 


1.045 



Inclusive jet cross section in bins of Q 2 



1 


113.73 


1.8 


6.9 


4.9 


4.8 


3.2 


0.6 


0.4 


3.2 


1.062 


2 


94.06 


1.9 


7.1 


5.1 


5.0 


3.5 


0.2 


0.3 


3.2 


1.060 


3 


94.30 


1.8 


6.7 


4.8 


4.7 


3.1 


0.1 


0.0 


3.2 


1.063 


4 


85.41 


1.9 


6.2 


4.5 


4.3 


2.6 


0.1 


0.4 


3.1 


1.067 


5 


87.81 


1.9 


6.5 


4.7 


4.5 


3.0 


0.2 


0.7 


2.9 


1.071 


6 


4.24 


8.2 


15.1 


12.1 


9.1 


7.0 


5.0 


1.1 


2.2 


1.061 



Table 1 : Results of the inclusive jet cross section measurement using the inclusive k T algorithm 
in the Breit frame for the phase space 0.2 < y < 0.7, 7 < E T < 50 GeV and -1.0 < 
?7 Lab < 2.5. The multiplicative hadronisation correction factor as applied to the NLO calculation 
is shown in the last column. The contribution of ±1.5% from the luminosity measurement 
uncertainty is included in the total correlated uncertainty. 
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Normalised inclusive jet cross section in bins of Q 2 and E T 



total total single contributions to correlated uncertainty 



bin 


normalised 


statistical 


total 


uncorrected 


correlated 


model dep. 


positron 


positron 


HFS hadr. 


hadronis. 




cross 


uncert. 


uncertainty 


uncertainty 


uncertainty 


detector corr. 


energy scale 


polar angle 


energy scale 


correct. 




section 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


factor 


1 a 


0.168 


2.1 


5.4 


4.0 


3.7 


2.8 


0.6 


0.3 


1.8 


1.076 


1 b 


0.074 


3.1 


6.4 


4.9 


4.1 


2.0 


0.8 


0.6 


3.1 


1.035 


1 c 


0.015 


7.0 


9.5 


8.3 


4.7 


1.0 


0.4 


0.4 


4.3 


1.032 


1 d 


0.002 


18.9 


20.5 


19.7 


5.6 


0.0 


1.1 


0.3 


5.3 


1.065 


2 a 


0.184 


2.2 


4.7 


3.5 


3.2 


2.2 


0.6 


0.2 


1.6 


1.075 


2 b 


0.092 


3.1 


7.9 


5.9 


5.3 


4.0 


0.8 


0.4 


2.9 


1.034 


2c 


0.020 


6.8 


9.5 


8.2 


4.8 


0.1 


1.2 


0.4 


4.4 


1.041 


2d 


0.003 


18.9 


21.3 


20.1 


7.0 


3.1 


0.8 


1.1 


5.9 


1.044 


3 a 


0.199 


2.2 


4.2 


3.2 


2.8 


1.9 


0.5 


0.1 


1.3 


1.085 


3 b 


0.110 


2.9 


7.1 


5.3 


4.7 


3.5 


0.5 


0.0 


2.7 


1.032 


3c 


0.028 


6.0 


9.2 


7.7 


5.1 


2.4 


0.3 


0.2 


4.2 


1.030 


3d 


0.003 


18.1 


20.7 


19.4 


7.0 


1.3 


0.9 


1.2 


6.5 


1.039 


4 a 


0.228 


2.3 


3.8 


3.0 


2.4 


1.3 


0.5 


0.1 


1.2 


1.093 


4 b 


0.126 


3.0 


6.7 


5.1 


4.3 


3.3 


0.6 


0.3 


2.3 


1.035 


4c 


0.040 


5.4 


9.3 


7.5 


5.5 


2.7 


0.8 


0.7 


4.4 


1.025 


4d 


0.008 


13.3 


15.5 


14.4 


5.7 


2.7 


0.1 


0.3 


4.8 


1.035 


5 a 


0.239 


2.4 


3.7 


3.0 


2.3 


0.3 


0.7 


0.2 


1.5 


1.103 


5 b 


0.168 


2.9 


4.5 


3.6 


2.6 


1.1 


0.7 


0.2 


1.6 


1.040 


5c 


0.066 


4.5 


11.4 


8.6 


7.5 


6.5 


0.6 


0.1 


3.3 


1.038 


5d 


0.015 


10.3 


14.8 


12.7 


7.5 


4.9 


0.4 


0.3 


5.5 


1.046 


6 a 


0.225 


10.8 


11.7 


11.2 


3.2 


0.4 


1.8 


0.6 


1.7 


1.083 


6 b 


0.154 


13.4 


20.6 


17.3 


11.2 


9.5 


4.7 


1.9 


1.8 


1.050 


6c 


0.095 


17.8 


25.2 


21.7 


12.7 


9.6 


7.6 


2.2 


2.2 


1.029 


6 d 


0.055 


26.8 


33.5 


30.3 


14.2 


13.2 


4.0 


0.8 


2.9 


1.029 



Table 2: Results of the normalised inclusive jet cross section measurement, i.e. the average 
number of jets with 7GeV < E T < 50GeV and -1.0 < r] Lah < 2.5 per NC DIS event for 
the phase space 0.2 < y < 0.7. The multiplicative hadronisation correction factor as applied to 
the NLO calculation is shown in the last column. The Q 2 , Et binning scheme is displayed in 
Table [B 



15 



Inclusive Jet Cross Section 
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Figure 2: The single differential cross section for inclusive jet production as a function of Q 2 
(left) and of E T (right). The data, presented with statistical errors (inner bars) and total errors 
(outer bars), are compared with the results of NLOJET++, corrected for hadronisation and Z° 
boson exchange. The bands show the theoretical uncertainty associated with the renormalisation 
and factorisation scales and the hadronisation correction. In addition to the differential cross 
section, the ratio R = cr data / a th cory is shown. The band around R = 1 displays the relative 
error of the theory calculation. 
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Figure 3: The double differential cross section for inclusive jet production as a function of E T 
for six regions of Q 2 . The data, presented with statistical errors (inner bars) and total errors 
(outer bars), are compared with the results of NLOJET++, corrected for hadronisation and Z° 
boson exchange. The bands show the theoretical uncertainty associated with the renormalisation 
and factorisation scales and the hadronisation correction. In addition to the differential cross 
section, the ratio R = crdata/^thcory is shown. The band around R = 1 displays the relative 
error of the theory calculation. 
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Figure 4: The normalised inclusive jet cross section as a function of Et for six regions of Q 2 . 
The data, presented with statistical errors (inner bars) and total errors (outer bars), are com- 
pared with the results of the NLOJET++ and DISENT programs, corrected for hadronisation 
effects and Z° boson exchange. The bands show the theoretical uncertainty associated with the 
renormalisation and factorisation scales and the hadronisation correction. In addition to the dif- 
ferential cross section, the ratio R = a da ta/ 0"theory is shown. The band around R = 1 displays 
the relative error of the theory calculation. 
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Figure 5: Results for the fitted values of a s (Er) using the inclusive jet cross section for six 
regions of Q 2 . The error bar denotes the uncorrelated experimental uncertainty for each fit- 
ted value. The solid line shows the two loop solution of the renormalisation group equation 
evolving the averaged a s (Mz) from all determinations, with the band denoting the correlated 
experimental uncertainty. 
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Figure 6: Results for the fitted values of a s (ET) for six regions of Q 2 using the normalised 
inclusive jet cross section. The error bar denotes the uncorrelated experimental uncertainty 
for each fitted value. The solid line shows the two loop solution of the renormalisation group 
equation evolving the averaged a s (Mz) from all determinations, with the band denoting the 
correlated experimental uncertainty. 
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oc s from Norm. Inclusive Jet Cross Section 




Figure 7: Results for the fitted values of a) a s (/i r = E T ) averaged over all Q 2 regions, and 
b) a s (fi r = Q) averaged over all E T regions. The error bars denote the total experimental 
uncertainty for each data point. The solid curve shows the result of evolving a s (M z ) averaged 
from all Q 2 and E T regions, with the band denoting the total experimental uncertainty. The 
world average from PDG is also shown. 
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